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Abstract 

Block Toeplitz and Hankel matrices arise in many aspects of appli- 
cations. In this paper, we will research the distributions of eigenvalues 
for some models and get the semicircle law. Firstly we will give trace 
formulae of block Toeplitz and Hankel matrix. Then we will prove 
that the almost sure limit 7^,™^ il^'^) of eigenvalue distributions of 
random block Toeplitz (Hankel) matrices exist and give the moments 
of the limit distributions where m is the order of the blocks. Then we 
will prove the existence of almost sure limit of eigenvalue distributions 
of random block Toeplitz and Hankel band matrices and give the mo- 
ments of the limit distributions. Finally we will prove that 7^™) (7^,™-*) 
converges weakly to the semicircle law as m — )■ 00. 
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1 Introduction 



In random matrix theory, a very important object is the eigenvalue distri- 
bution of a random matrix. If ^ = {aij{uj))f-^-^ is a real symmetric random 
matrix where the (ajj(a;))'s are random variables on a probability space VL 
with a probability measure P, then the eigenvalue distribution of A is 
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where Aj(a;)'s are the N real eigenvalues of A. 

The asymptotic behavior of the eigenvalue distribution is of mueh im- 
portance. In [23,24], Wigner got the semicircle law for a wide class of real 
symmetric random matrices and this great result caused much development 
of random matrix theory. Recently in a review paper (see [1]), Bai pro- 
posed the study of random matrix models with certain additional linear 
structure. The properties of the distributions of eigenvalues for random 
Hankel and Toeplitz matrices with independent entries are listed among the 
unsolved random matrix problems posed in [1]. In [5], Bryc, Dembo and 
Jiang proved the existence of limit distribution 7^ and 7^ of real symmetric 
Toeplitz and Hankel matrices. The moments of 7^ and 7^ are the sum of 
volumes of solids. Hammond and Miller in [9] also proved the existence of 
7j, and 7^^ independently. In [15], Liu and Wang proved the existence of 
limit distribution 7^ and for real symmetric, complex Hermitian band 
Toeplitz and real symmetric band Hankel matrices. We notice that Basak 
and Bose (see [2]) and Kargin (see [13]) also did the same work indepen- 
dently. Especially, the limit distribution of random Toeplitz band matrices 
with bandwidth = o{N) is Gaussian. 

Block Toeplitz and Hankel matrices arise in many aspects of mathemat- 
ics, physics and technology (see [8,10,14,19]). A block Toeplitz matrix is a 
block matrix which can be written as 
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where Ag = {aij{s))'^j^i is an m x m matrix, Vs G {—N + 1,...,N — 1}. 
In [7], Gazzah, Regalia and Delmas researched the asymptotic behavior of 
the eigenvalue distribution for block Toeplitz matrices. In [18], Rashidi Far, 
Oraby, Bryc and Speicher have proved the existence of the limit distribu- 
tion of eigenvalues for a random block Toeplitz matrix whose blocks are 
selfadjoint m x m matrices as m — t- 00, which implies that the double limit 
lim lim gives the semicircle law. In this paper, we will study the limit 

distributions of real symmetric random block Toeplitz and Hankel matrices 
as lim and lim lim . 

In a block Toeplitz matrix with the form mentioned above, we suppose 
the ajj(s)'s are real random variables. For symmetry, we need As = {A-g)^ . 
For independence of the elements, we suppose: 
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(1) aijj^(si) and 042^2(^2) are independent if 7^ IS2I5 

(2) If s / and (hji) / («2,i2) then ai-j-^{s) and independent, 

(3) If / (i2,j2) and (ii, ji) ^ {32,12) then aiiji(O) and 012^2(0) are 
independent. 

In addition, we need the following uniform boundedness conditoin: 
(4) 

E{a^j{s)) = 0, E{\ai^{s)\^) = 1,-{N - I) < s < N - \,1 < i,j < m (1.1) 
and 

{|aij(s)|^|l <i,j<m} = Ck,m < +00 (1.2) 



sup 

jv e N 

-{N -l)<s<N-l 



and 



sup Ck,m = Ck < +00. 
meN 

A block Hankel matrix is a block matrix which can be written as 



(1.3) 



( ^AT-l A.N-2 

^N-2 ^7V-4 
^W-S 



V ^0 



A. 



A, 



^0 

^-2 

^-(JV-l) / 



where = (ajj(s))™^;^ is an m x m matrix, Vs G { — iV + 1,...,A^ — 1}. 
Similar to random block Toeplitz matrices, we suppose the aij(s)'s are real 
random variables, Ag = (^s)^- In addition, we assume: 

(1) aj-^jj(si) and ai^j,2^{s2) arc independent if s\ 7^ S2, 

(2) If (ii,ii) / (i2,j2) and (ii,ji) / (^2,^2) then ai^j^is) and aj2j2(s) are 
independent, 

(3) 



E{a,j{s)) = 0, E{\aij{s)\^) = 1, -(iV - 1) < s < iV - 1, 1 < i,j < 



m 



and 



and 



sup 

jv e N 

-{Af -l)<s<A'-l 



{|ajj(s)|''|l < i,j <m} = Ck,m < +00 



sup Ck,m = Ck< +00. 
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We will firstly give trace formulae of block Toeplitz and Hankel matrices 
in Section 2. Using those trace formulae, we will prove that the almost sure 
limit distributions of block Toeplitz and Hankcl matrices exist in Section 3. 
In that section we will also give the moments of the limit distributions. In 
Section 4, we will prove the existence of almost sure limit distributions of 
block Toeplitz and Hankel band matrices and give the moments of the limit 
distributions. In Section 5 we will prove that for block Toeplitz and Hankel 
matrices, the double limit lim lim gives the semicircle law. 

m->-oo N-^oo 



2 Trace Formulae of Block Toeplitz and Hankel 
Matrices 

Definition 2.1. Let T be an mN x mN matrix and consist of N'^ blocks. 

If T has the form 
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where {A_ 



(iV-l)) ^0, 



,An-i} is a set of m X m matrices and Aij 



Ai-j = {apq{i — j))^g=i , then we call T a block Toeplitz matrix. 

Let H be an mN x mN matrix and can be written as H = ^>T where T 
is a block Toeplitz matrix and 



$ = 



/ 








Im \ 





0/ 



where Im is the mx m unit matrix. Then we call H a block Hankel matrix. 
For convenience, let In = N — 1. 

Lemma 2.2. Let T = (^i-j)f^=i be a block Toeplitz matrix and Ag = 
(apq(s))^g^^ where —bN < s < bN', ^ < P,q < th. Then we have a trace 
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formula 



N bN k I 

i=l h,-,jk=-bN 1=1 «=1 

N bjsi m k I 

= E E E (^titAh)---at^ti{jk)Yll[i,N]{i + J2^i^^. * 

j=l ji,...jfc=-6Ar ti,---,tfc=l ^=1 9=1 



0, E j, 

9 = 1 



Lemma 2.3. Let H = $T 6e a block Hankel matrix and T = {Ai^j)fj^-^; As = 
{0'pq{s))^^g=i where —bN < s < bN] 1 < p,q < m. Then we have the trace 
formula 

tr(ij'=) = 



( N 


bN 
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E 




•Jfc= 


' N 
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»=1 jlr 





in 

1=1 
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E( 

9=1 

E( 

9=1 



0, E {-i)'^J<; 

9=1 



5 fe odd 

2i-l-N, E (-1)«J9 

9=1 



where 



tr(^ii • • • J = E (ii) • • • «tfeti Uk)- 

ti,...,tk=i 

To prove the above lemmas, we consider Kronecker product of two ma- 
trices (sec [12]). Let A = {a,ij) be an m x n matrix and -B be a p x g matrix. 
The Kronecker product of A and B is an mp x nq matrix: 



/ aiiB ■ ■ ■ ainB \ 



B 



\ Oml-B • • • flmn-B / 

Proof of Lemma 2.2. Let B and F be two N x N matrices and 
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1 
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0/ 



- {^i,j+l)f,j=l 
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1 0/ 
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Then T = ® A_j + Y,F^ ®Aj. Let e\ be the ith unit vector in 

and 62 be the jth unit vector in M"*, then we have 

(g) A_i{e{ (g) e^) = {B^e[) ® (A.^e^) = /[i,7v](i - 04^' ® 

and 

thus 

i=0 /=1 



so 

r 

T\e\ ® e^2) = E n hm + E ^9)^1 ® (^'^ • • • ^h4)- 

h,-.-,lk=—bN '''=1 q=l 

As {e| (g) 62 I 1 < i < AT; 1 < j < m} is a standard basis of W^^, we have 
the trace formula 



tr(r^') 




A' m 






i=i j=i 
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1=1 h,.. 


■I'fc— 


N 
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1=1 h,.. 


•>'fc= 


N 
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1=1 h,.. 


■,h= 



?=1 j=l 



m 



k 



q=l '^fi ' J = l 
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For tT{Ai^ ...Ai^)= tr(Aii ■ ■ ■ ^/J, we get 

N bN k I 

i=l n,...,jk=-bN 1=1 q=l 



Note that 



ti,---,tk=i 

and we directly get 

N bN m k I 

tT{T'') = J2 J2 ""tit, (ii) • ■ • at,t^ Uk) n ^[i,m 31^^ ^ . ■ 

1=^ ji,-,jk=-bNti,...,tk=l 1=1 q=l ^'^1^" 

□ 

Proof of Lemma 2.3. ^ = P ® where Im is the m x m unit matrix and 
P is an N X N matrix and 



P = 



( 


V 1 



1\ 

1 

y 



So if A is an X matrix and i? is an m x m matrix, then we have 
^{A ®B) = {P® Im){A (^B) = (PA) tE) B. Note that Pe\ = ef +^-\ 
As in the proof of Lemma 2.2, we have 

me\<S>ei) = iP0lm)i Yl I[i,N]{i+l)e\^'<S>{A'ei))= Yl %,iv](^+/)ef 

l=—bN l=—bN 



and then 



bff k r 

E n^[i,iv](^-E(-i)%K 

ji,-,jk=-bN r=l 9=1 



i- E (-1)''J<; 

9=1 



^ {Ai^^ ■ ■ ■ Ai^)e'2 k even, 



6jv k r 

E n%,Ar](^-E(-i)%K 

,ji,-Jfc=-6Ar '■=1 9=1 



N+i-{i- E 



(^^^•••A^Je^ A; odd. 
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= < 



Thus we get 

tr(i?'^) 
N m 

i=i j=i 

' N bM k I 

E E tr(A,,.--^,-jn%,Ar](^- E(-iri.)'5 . ^even 
i=ln,...,jk=-bN 1=1 9=1 o,j:i-i)ijg 

9=1 

iV bjv fc « 

E E tr(^.V^.-jn^[i,iV](^- E(-l)^J<z)'^ . ^odd 

9=1 

where 

m 

ti,---,*fc=i 

□ 

3 Limit Distributions of Random Block Toeplitz 
and Hankel Matrices 

Let [ty denote the number of elements in an arbitrary finite set V. We now 
review the concept of partition (see [15]). Let [n] = {1,2, ...,n}. 

(1) We call TT = {Vi,V2, a partition of [n] if (j Vj = [n] and Vi f] Vj = 

3=1 

if i ^ j. 

(2) For Vp G [n], we define 7r(j)) = i \ip G Fj. We write p 9 if 7r(p) = '^('z)- 

(3) Let V{n) denote the set of all the partition of [n]. We define V2{n) = 
{vr = {Vi, V2, Vr} G r{n) \ Wi = 2, Vi} and Vl{n) = (tt = {F, V2, K} e 
7^2 (?^) I Vi has exactly one even element and one odd element, Vi}. 

Definition 3.1. Let G N, tt = {Vi, V2, Vk] G ^2(2^;), Vr = {a^, 6^} (1 < 
r <k), then tt determines a system of linear equations with unknown vari- 
ables ii,i2, (set t2k+\ =ti): 

Ua, = tb.+i l<i<k ^^^^ 
\th^ = l<i<k 

The number of linearly independent equations of this system is denoted by 
/(tt). Clearly we have < /(tt) < 2k. 
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Theorem 3.2. LetT^q = {Ai-j)f^j^-^ be anmNxmN random block ToepUtz 
matrix where As = (ajj(s))™=i, ^-s = {Ag)'^ , andaij{s) 's are random vari- 
ables satisfying the four conditions mentioned in introduction. Let = 
-^=^, then eigenvalue distribution fi^^ converges almost surely to a sym- 
metric probability distribution jj^^ which is determined by its even moments 



„ 2k j k 



where 



I 1 if q is the smaller element of K-c„i, 
eAq) = { ^' (3-2) 

1—1 if q is the larger element of ^-(5) • 

Proof of Theorem 3.2. The fcth moment of jj,^^ is given by 



Using Lemma 2.2, we get 

^ N bM m k I 

^^^^ ^ ^ ^ -£^(atit2(ji)---atfcti(jfc))n%-^](^+IZ-^«)'^ " . 

i=lji,...,jk=-bNti,-,tk=l 1=1 q=l ^^i^" 



We will calculate the limit lim mk n- 
Observe that 

\E{atit2{ji) ■ ■ ■ atkti{jk))\ <Dk (Vji, jfc, ti, t^), 

where = (max{Ci, Ck})''. If \jp\ / \ jq\, then at^tp+i Up) and at^tg+i Uq) 
are independent. Let p G {1, 2, k}. By the independence conditions and 
E{aij{s)) = 0, we observe that if E{at^t2{h) ■ ■ ■ atktiUk)) 0, then for any 
p G {1, k}, 3q G {1, A;}, q ^ p, s.t. \jq\ = \ jp\. So there are at most [|] 
different elements in y = Further, we have 

^(atit2 (ji) • • • «tfcii (ifc)) = I XI XI ^("*i*2 (Ji) • • • otfeii Uk)) 
ji,—,jk=-bN s=liV=s 

^ ^^^'^(6^ + 1-.)!'^'= 
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where = ^{tt = {Ui,...,U,} G V{k) | > 2,Vi}. 
So 



6jv 



E E{at,tM---atktAjk)) =0(iV[t]) and then mfe,jv = 0(Ar[|l-t). 

Jl>---Jfc=-&JV 

Thus for odd fc, m^^Ar = o(l). Now we only have to consider m2k,N- 

TT is a partition of [2fc] = {1, 2, 2k} and p ~,r 9 |jp| = then 

we have 



UM+l iP[i-JV](^ + 2^.7?)^ 2. 

j=l7re7'(2fe) n.-.32k=->>N ti,...,t2k=l ^ ' l=\ ?=1 'i'l-''' 

+ o(l). (3.3) 

Again by the assumptions in the introduction, the contribution of the par- 
titions which are not pair partition to m2k,N is o(l). So we only have to 
consider the pair partitions. Suppose tt G V2{2k). 11 p q, then jp = jg or 

2k 

jp = —jq. Under the condition ig = according to (3.3), considering the 

q=l 

main contribution to the trace, we should take jp = —jq. Otherwise there 
exists po,qo € [2k] such that 



^ 2k 

jpo — jqo — 2^jpo + jgo ~ Jg)- 



9=1 

We can choose other k — 1 distinct numbers, which determine jp^ = jg^ and 
then there is a loss of at least one degree of freedom and the contribution of 
such terms is 0{N~^). Therefore we have 

'm2k,N 

1^ 1^ 1. (mN)k+^ ii%,iv](^ + Z^J9) 

i=l 7re7'2(2fc) Jl.- -.J2fe=-''JV ti,...,t2fe=l ^ ' 1=1 q=l 

+ o{l) 

Ev^ ^{O'titi [Jl) • • • at2kti{j2k)) jj T A \ 

1^ 1^ 1^ MvF+i ivhmy' + z^h) 

i=l 7re'P2(2fe) 31, ■■■•hk = -bN tl,...,t2fc=l 1=1 ?=1 

h ^ 0, vt 

+ oil). 
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Now we have to calculate 

m 

XI ^(atit2(ii)---at2fcti(i2fc)) 

with the assumption that jt 7^ 0, Vf. For convenience we let t2fe+i = ii- 
Suppose 

then for any p ~^ g we have 0*^4^+1 (jp) = at^tq+Aiq)- Suppose p < q. 

• O'tptp+l Up) ■ ■ ■ (Hqtq+l (jq) ■ ■ ■ 0,t2ktl (i2fe)) 

= 1 • E{at^t2 (ii) • • • Op) • • • %t,+i (i?) • • • «t2fcti ihk)) 

where ^ marks omitted index. Note that jp = —jq, Aj^ = (Aj^)'^ and Aj^ 7^ 
Aj^ because jp = —jq ^ 0. at^tp+i (jp) lies in the tp-ih. row and the fp+i-th 
column of Aj^ and atqtq+iijq) lies in the tg-th row and the tg+i-th column of 
Aj^. So tp = tq+i, tq = tp+i. As TT = {Vi, V2, Vr} where Vi = {a^, 6^} (1 < 
i <r), we have a system of equations: 



ai 



tbi+l 
tai+l = hi 

= h. 



By Definition 3.1, there are /(vr) linearly independent equations in this 
system. So there are 2k — /(n) variables taking values freely in {1, 2, m} 
and the number of solutions of this system of equations is rn'^''~f^'^^ . In other 
words, there are m'^'^~^^'^^ different (ti, ...,t2jk)'s such that 



and this implies 
So 



E{atit2{ji) ■ ■ ■ at2kti{hk)) 7^0 
E{atit2Ui) ■ ■ ■ at2kti{j2k)) = 1- 



*lr--)*2fe = l 
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Thus we have 



m2k,N 

^ N bN 2k I 

= T^^E E E n^M('+E^.) 

^ i=l 7re7'2(2fc) ii,. ■•,i2fc = -i>iv 1=1 9=1 

n 7^ 0, vt 

P ~ir 9 Jp = -jg 

+ o{l) 

N bN 2k I 

= (SD^E E E n^'^t'+E^.) 

^ ' i=l 7re7'2(2fc) Ji. - J2fc=-''jv 1=1 9=1 

P~7r9 -!=^ jp = —jq 



Now for any r G {l,...,fc}, let = Ja^, then 



,Tr ii q = a. 
—Xr if g = 6, 



Remember Eq.(3.2), and then we have jq = i '^^^^ " '^^^^ , so 

(^-x^(g) ife^(g) = -l 

J? = ^n{q)x.„{q) and 

iV bjv 2k I 

m2k,N = {mN)-^-'Y, E E "^''"^^'^^n%.iV](^+E^-(^)^-w)+'^(l)' 

i=l 7reP2(2fc) a;i,...,Xfc=-fejv '=1 9=1 

thus 

» 2fc i k 

N-^oo ^ — ' JiQ iiv[-l life — ^ 



Letm2fc= E Y\l[Q,i]{xo + y2^-^^^^^^ii)'>T\'^^'- 

.eP2(2fc) •^[o,i]x[-i,i]'=^=i g=l 



and m2fc-i = 0(A; G N), then for any A; G N we have hm n = fnk- It 

N—^oo ' 

is easy to see that m2k < X] w}'^^~^^'"^ < {2k — 1)!! • m''~^ and then 

7reP2(2fc) 

using Carleman's theorem (sec [6]) wc know that the hmit distribution 7^"*-* 
is uniquely determined by its moments {mk}^^Q- 

Now we prove the almost sure convergence. It is sufficient to prove 

00 -. 

J2 J^E[{tvX'^ - E{tvX'^))') < 00. (3.4) 

N=l 
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trX^ = {mN)-2tiTj 



N 



N 



k 



= {mN) (^tit2ih)---atktAjk)Yll[i,N]ii + Yjq)^ k 

0, E jq 
9=1 



1=1 h,...,t);=l ji,...,ji.=-bN 



1=1 



q=l 



For convenience we let 

3 = {ji,-,jk), t = {ti,-,tk), 

k I 

y4[i,t,j] =otit2(ii)---atfcti(jfc) JJl[i,iv](i + ^jg)(5 k , 

1=1 q=l 

N m biq 

and we use ^ to denote S j then we have 

i,t,3 i=l tl,...,tli: = l ji,...,jk=—bN 

1 „//. Tpfi.„vk 



l^E(^{mN)--2 A[i, t, j] - E{{mN)-'2 ^ A[i, t, j]))^ 



= ]v^^^(E E ^(^[^^t^j1-w,t^j1)) 

N m bff 

where E = E E E • 

For given = (jj, f = = Ut-Jllf = Uh-Jt), 

let 

J = Ul, ...,jf, ...Jt) e {-bN, i'Ar}"^ 

Set S"^ = I, Ij^I, Iji I) bfcl}- We use p{3) to denote the number 
of different elements of S'^ . 

We construct aset of numbers with multiplicities Sj = {jl, •••,J^, ...,J^}. 
Please note that if = = a, then a appears twice in Sj. Let 

Si,..., 5'p(j) be subsets of Sj such that: 

(a) for all w, the elements of Su, have the same absolute value; 

(b) if wi ^ W2, then the absolute value of the elements in Sy,i is different 
from the absolute value of the elements in ! 
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p 

(c) IJ Sr = Sj. It is easy to see that ^i, S'p(j) are uniquely determined 

r=l 

by J. Then we have 

{-bN, 6jv}^'' = WiUW2UW3 

where 

Wi = {3e {-bN, bN}^''\p{3) <2k- 2}, 

W2 = {3e {-bN, bNj'^^lpiJ) = 2k- 1}, 

W3 = {Je {-bN, bN}'"'\piJ) = 2k}. 
Therefore we have 

^s((trX^-£;(trX^))') 



I ,1 ,1 ,1 



.1 .2 ,3 ,4 JgVl/^ v=l 
4 



+ ]^^2^^( E E ^(^[^^t^J1-w,t^j1)) 



tl,t2,t3,t* 



JeVF2 1^=1 



Jew3v=i 

t^t^t^t* 

= $i + $2 + ^3- 

For given if,..., j4,ii,...,i4^ii,...,i4^g^ppoge^(j^(^[^.^t.j.j_^^^^^ 

v=l 

k 

0. Then = (1 < v < 4) and from independence conditions we know 

i=l 

that: 



\/j:,3{uuvi)j^{u,v)s.t. \j:\ = \j:i\ (3.5) 

and that 

Vui, 3i;2 7^ VI and 'Ui,n2 s.t. \ = (3.6) 
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otherwise, 



E{ll{A\i\es]~EA\i\es])) 

v=l 

= E{iA[i-\t'^\r]-EA[i-\t-\r]))-E{ n iA[i\eS]-EA[i\t\r])) 



1 <v < 4 

V ^ VI 



We now evaluate ^>2 and ^>3. 

Evaluation of $i. Suppose p e N and p < 2k — 2. There are 
at most Rk,p ■ {b^ + 1)^ • 2^^^ different J's that satisfy p{3) = p, where 
Rk,p = tt{^ = {^1' ••' ^pik ^ ^W}- So we have 

p=l 



where 



4 

M = sup I max l^^fTT t^ j''] - £;A[^^t^^])]|| 

Je {-6„,...,6jv}'"' ^^=1 
t" S {1, ...,m}'= 

and C is independent of A^. From (1.1), (1-2) and (1.3), we know that such 
M exists. 

Evaluation of $2. Suppose J G W2. For any v G {1,2,3,4}, if 
3u G {1,2,..., k} such that Vui / u, / linj) then and its absolute 

A: 

value are determined by jl, ...Jl E J'i' = 0- Thus 

1=1 

A: 

tt{j = (ii\...,jt)eW^2|Vi;,5^ir = 0and ^i^J s.t. Vus 7^ m, b^J 7^ 

< Qk-{bN + lf'-^-2^K (3.7) 

where Qfc = tt{vr = {^1, V^2fc-i} |vr G 7^(4A;)}. When J G W2, there are 
two situations. One is that 3ii,i2 such that ^Si-^ = '^Si^ = 3 and for any 
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i ^ {^l,^2}, tt'S'i = 2. The other situation is that 3zi such that ftSj^ = 4 and 
for any i / ii, ftSj = 2. In the first situation, suppose Si^ = {jxdydz}^ 
~ {jw^Ju^jv}- If a = b = c and d = e = f, then wc can find g G 
{1, 2, 3, 4}\{o, d}. Then for any r, there is one unique element which has 
the same absolute value as jr in Sj. Then from (3.6) we know that in 
{jf, there is at least one element which has a different absolute value 

from the others and then (3.7) is obtained. If a ^ b and a ^ c, then has a 
different absolute value from the other elements in {jf, and then (3.7) 

is also obtained. In the second situation, suppose Si^ = {j^jjy^j^jt}- 1^ 
{a, 6, c, d} 7^ {1, 2, 3, 4} then similarly to the first situation we can get (3.7). 
If {a, 6, c, d} = {1, 2, 3, 4}, then jf has different absolute value from the other 
elements in {jf, J^i thus wc also get (3.7). From the above discussion we 

4 

know that for any J G such that E{ \[ (A[i^ t^, j'^] - EAIV', t^, j'^])) ^ 0, 

v=l 

(3.7) is satisfied. Then we get 
$2 

3&W2 v=i 

ti,t2,t3,t* 

< ]v^^^'-''Q^(^iV + l)^^-^2^^M 

< D-b-^^ 

where D is independent of A'^. 

Evaluation of $3. When J G W3, we have (tS'i = 2 for all 

i G {\, ...,2k}. Like the situations mentioned above, we have: 3ui G 
{l,2,...,fe} such that \lu2 7^ ui, |j^J 7^ |ju2l' then j^^ and its absolute 

k 

value are determined by if, for j| = 0- Suppose 



1=1 

\fu^ = and / G {l,2,3,4}\{l,e}, then 3u4 G {1,2,..., A;} such that 
7^ -^4, Ij^gl 7^ iJu^l, then and its absolute value are determined by 

k 

ji, -Jli-ijli+i, -Jk foi' E Ji = 0- For jl^ and j4 are determined, we 
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have 

4 

tt{j = e W^\E{l[{A\i\t\r]-EA[i-,e,r])) + 0} 



< (6jv + l)^~'-2^'=-r2(4fc) 



and thus 
$3 



= j^2^A J2 E^(^(^[^^t^j1-w,t^r]))) 

JeWs v=i 



t^t^t^t* 



where F is independent of N . Finally, 

^E{{tvX% - E{\,TX%)f) = $1 + $2 + $3 < (C + D + F) • 6-2. 
For In = N — 1, (3.4) is proved and we have proved the theorem. □ 

Remark 3.3 (Hankel block matrices and complex Toeplitz case). If 

H]\f = (^7v+i-i-j) is a block Hankel matrix which is defined as in the in- 
troduction, then the eigenvalue distribution of Hn/VitiN converges almost 
surely to a distribution 7^*^^ which is determined by its even moments 

2k j k 

".-(7ir')= E w / , , ,„.nvi(-c-i:(-i)'-.(,))n<i-. 

7r6Pl(2fe) ^10,1JX[-1,1J j=i q=i 

(3.8) 

where r(m,7r) = tj{(ti, .., t2fe) | E{at^t2{ji), ■■■,at2kti{j2k)) O}. 
Suppose TT = {{ai,bi},...,{ak,bk}}, then 

r{m, tt) 

= tt{(*i' ■•'*2fc) |1 < ta„h, < m; 1*"^ = or |*«^ = (1 < z < fc)} 

= tt{(*l' - -'^Sfc) I 1 < iai,i&i < "i; Ui = *(6i+l); *(ai+l) = (1 < ^ < ^)} 

+ 0{m'') 

= m^'^-^W + 0(m'=) (3.9) 
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If the blocks of i^jv are Hermitian matrices, the results would be similar. If 
the blocks in a block Toeplitz matrix T = are complex matrices 

and A-s = i^s)'^, then the results would also be similar. 
Remark 3.4 (relation with dynamical system), /(vr) has a relation with 
dynamical system. Let tt = {Vi, V2, V/j}, = {ar!^r}(l < f < k), 
then /(vr) is the number of linearly independent equations of (3.1). Now we 
consider a discrete dynamical system. Give [2fe] = {1, ...,2A;} the discrete 
topology. Let (/> be a self-homeomorphism on [2/c] and (j){ai) = hi + l, (/>(6j) = 
+ 1 (1 < i < k:). Consider a continuous map : Z x [2k] — )■ [2k] such 
that ip{s,t) = 4>''^{t) {s e Z, t e [2k]) then ■0 becomes a dynamical system. 
Obviously we see that the number of orbits of this dynamical system equals 
the number of independent variables in (3.1) and then equals 2k — /(tt). 

4 Block Toeplitz and Hankel Band Matrices 

Definition 4.1. Let = (Ai_j)^-=i be an mN x mN block Toeplitz 
matrix. We call T a block Toeplitz band matrix if 3 67V < s.t. As = 
when |s| > 6jv. We call bN the bandwidth of the matrix. 

Theorem 4.2 (proportional growth). LetT^ he a block Toeplitz band matrix 
with the bandwidth bjy ^ bN, b G (0, 1]. Take the normalization 

Xn = TN/^m{2-b)bN. 

With the notation and assumptions of Theorem 3.2, //^^ converges almost 
surely to a symmetric probability distribution 7^"*^ (6) which is determined 
by its even moments 



Ik j k 

m2fe(7f)(&)) = 7^^ E rn'-'-^(-^ [ lll[o,i]{xo+b^e^{i) x^^^)l[d 

where e-nii) = lifiis the smaller number o/K-(i) ^^'^ ^wii) = ~1 otherwise. 

Theorem 4.3 (slow growth). Let Tjv be a block Toeplitz band matrices with 
the bandwidth bjy = o{N) but bjy ^ 00. Take the normalization 

XN = TN/{V^mbN). 
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With the notation and assumptions of Theorem 3.2, fi^^ converges weakly 
to a distribution 7^"*^ which is determined by its even moments 



7rG'P2(2fc) 

In addition, if there exist positive constants eo and C such that 

bN >C-N^+'°, 
then /i^^ converges almost surely to 7^"*) . 

Proof of Theorem 4-2. It is easy to see that Lemma 2.2 is also right for block 
Toeplitz band matrices but now 6jv is no longer N — 1. Let 



mk,N = I x'di,,^ = :^E{i^X%) = ^(m(2 - b)bN)-'i E{iTT, 
Using Lemma 2.2, we get 
rrikN 



1 E E ^'"■■■'';'>-'>-.fe» n^i..~i('+E^»)^, 

i.lj, y.=-J„(„...,(..l ("'AT). ,.1 I 

■^(j)'E E E ^'"■'''^' ■■,°'".fa» n%..|(.+E^-.)^ 



Similarly as in the proof of Theorem 3.2, we know that when k is odd, 

N 6jv c^/ / ■ ^ / ■ \\ ^ ' 

»=iii,...jfe=-6jvti,...,t.=i (miV)2+ o,E^j, 

and then m^^jv = o(l). 
For m2fc,Ar, 

AT bjv m „f , . s 2k I 

f^ATl-t+i [[^lhN]{^ + 2^h)S 21, 

i=l h,-,hk=-bNti,...,t2k=l ^ ' 1=1 9=1 \W 

N bis, 2k I 

1=1 TTe'P2(2k) Xl,...,Xk = -bN 1=1 9=1 
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So 

'm2k,N 

^ ^ f. N bN 2k I 

i=l TveV2{2k) xi,...,Xk=-bN 1=1 9=1 

thus 

\uam2k,N=-7^—^ E "^'"'"^^"w^^,, ^ nv](^o+6E«-(«)^-w)nd^'- 

^ ■KeV2{2k) ^[0,11x1-1,1] -^^ 

Let = lini ^2^. and m2A;-i = 0(/c G N), then for any G N we have 
hm mfe jv = JTife- It is easy to see that 

m2k < E ^^'"^"-^^''^ < (2^^ - 1)!! • m''-'^ 

7reP2(2fe) 

and then using Carleman's theorem (see [6]) we know that the hmit distri- 
bution ^^\b) is determined by its even moments {m2k}- Similar to the 
proof of Theorem 3.2 we can easily prove that /x-^^ converges almost surely 

to 7^"^) (6). Then Theorem 4.2 is proved. □ 
Proof of Theorem 4-3. Let 

mk,N 

= J ^'d/.,^ = = ^(2m5^)-tE(trT^) 



Similarly as in the proof of Theorem 3.2, we know that when k is odd, 
mk,N = o{l). For m2k,N, 

f^ATl-t+i [[^lhN]{^ + 2^h)S 2. 

i=l ji,...,j2k=-bNti,...,t2k=l ^ ' 1=1 q=l \W 

N 2k I 

= N-'^-'Y. E E ^'"'"^^"^n^[i-ivi(^+E^-(9Kw)+«(i)- 

1=1 TTe'P2(2k) Xl,...,Xk = -bN 1=1 1=1 
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So 

rrikM 



AT 6jv 2k I 



&JV (=1 9=1 



^ ' 1=1 Trer2{2k) xi,...,Xk=-' 

N bM 2k I 

jTWTw^ ^ ^ -'^-^-^^'^)n^[i>^i(^+E^^(«K(.))+«(i) 

i=i TreV2{2k)xi,...,Xk=-bN 1=1 ?=1 



Thus 

\xi 



lim m2,,iv = 4 E "^'"'"-^^"^ / fld.^ 
Ar-*-oo 2'^ ^-^ /[nilv[ 111* 



7reP2(2fc) 

Let m2k = lim n and m2k-i = 0(A; G N), then for any fc G N we have 
hm mfe jv = mfc. It is easy to see that 

N-KX> ' 

m2k < Yl rn'^-^-f^''^ < {2k - 1)!! • m'^"^ 

7r6P2(2fe) 

and then using Carleman's theorem (see [6]) we know that /x^^ converges 
weakly to a distribution -y^™) which is determined by its even moments 
{m2k}- Similarly as in the the proof of Theorem 3.2 we know that 



^s((trX^-^(trX^))') <B.6 



2 

TV 



where B is a constant and is independent of N. For bjsf > C ■ iV2+^o, we 
have 

^^((trX^ - E{tvX'^))') <B.C-'. 7V-i-2^o. 



Therefore 



oo 

j;^£;((trX^-£;(trX^))^) <oo 

7V=1 

and then //^ converges almost surely to the limit distribution j^^^ □ 
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Remark 4.4. In [15], the authors proved that the hmit of eigenvalue distri- 
bution for band Toephtz random matrix with bandwidth = o{N) is the 
standard normal distribution N(0,1). But in this paper, the matrix is block 
Toeplitz and the conclusion will be different. The expectation of 

^(™) is 

and the variance of 7^"*) is ^ m~^^'^^ = 1. The forth moment of 7^"*) is 

7reP2(2) 

^ m^~-^^'^^ = 2 + ^ which is not the forth moment of N(0,1) if m 7^ 1. 

7re'P2(4) 

So 7^"') is not N(0,1) when m 7^ 1. 

5 Convergence to Semicircle Law 

Suppose TT = {{ai, 61}, {ofc, 6^}} G V2{^k), the system of linear equations 
determined by tt is (see Definition 3.1) 

ta, = h,+i l<i<k 
hi = ta,+i l<i<k 

where t2k+i = ii- This system can be rewritten as 

*si(l) = *si(2) = • • • = tsi(ri) = isi(l) 

< . . . 

Jsp{l) = *sp(2) = • • • = isj,(r-p) = tsp{l) 

such that {si(l), si(ri)} U ••• U {sp{l), Sp{rp)} = {1, 2, 2A;} and 
{si{l),...,Si{ri)} n {sj{l),...,Sj{rj)} = $]ii ^ j. We call ts.(i) = ts^{2) = 
■ ■ ■ = tg.^^.-j = is.(i) a circle of vr. We use ^(Tr) to denote the number of the 
circles of vr. For' example, if tt = {{1, 2}, {3, 4}, {5, 6}, {7, 8}} G ^2(8), then 
the system of linear equations determined by tt is 



h 


= ^3 




= t2 




= t5 


U 


= U 




= t7 




= t6 




= h 




= t8 
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This system can be rewritten as 





= t3 






< t4 


= U 






M 


= t8 



So this TT has five circles: h = ts = = ty = ti, t2 = t2, ^4 = ^4, tQ = tQ and 
ts = ts, thus g{'K) = 5. These circles can be denoted by 







I t4 



For /(vr) denotes the number of independent equations of the system, it 
is easy to sec that g{7r) = 2k — /(vr) for any vr G P2(2fc). 

Lemma 5.1. For any tt G V2{2k), g{Tr) < A; + 1. 

The proof of Lemma 5.1 can be found in [22]. 

We now review the concept of noncrossing partition (see [21]). A parti- 
tion TT G V{n) is called noncrossing if whenever four elements 1 < a < 6 < 
c < d < n are such that a ~,r c and b ~,r d, then a ~,r b ~,r c d. 

Lemma 5.2. For any tt G 7^2 (2A;), g{7r) = k + 1 if and only if n is non- 
crossing. 

The proof of Lemma 5.2 can be found in [20]. 
Theorem 5.3. Suppose defined as in Theorem 3.2. As m — >■ oo, 



converges weakly to the semicircle law w{x), i.e., 



w{x) 



2n 





\x\ < 2, 
\x\ > 2. 



Proof of Theorem 5.3. From Theorem 3.2 we know that the odd moments 
of are all zero and its even moments are 

'T 

2k j k 

m2fe(7f))= J2 rn'-'-^^-'^ [ H ^[o,i](^o+^ 6,(g) Hd 

7reP2(2fe) [0,1] X [-1,1]%=! q=l i=o 



XI 
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where 

I 1 if g is the smahcr element of KrCo-) , 
1—1 if g is the larger element of T^^^) . 

From Lemma 5.1 and Lemma 5.2 we know that if vr G V2{'^k), then 
/(tt) = A; — 1 if TT is noncrossing, otherwise /(tt) > — 1. So we have 

„ 2fc i A; 

vrePaCsfc) •^[".i]xl-i'i] j=i q=\ ;=o 



TT IS noncrossing 



Prom [4] we know that when tt is noncrossing, 

» 2fc j k 

\ n-^[o,i](^o + y'e7r(g)a;7r(9))n^^^ = ^• 



So 



limm2fc(7f))= ^ L 



7reP2(2fc) 
TT is noncrossing 



Prom [11] we know that 

tl{7r G P2(2/i;)|7r is noncrossing} = Cfe 
where is Catalan number. So 

lim m2fe(7f)) = Cfc. 

Por Cfe is the 2A;-th moment of w{x) whose odd moments are all zero, we 
know that 7^"*^ converges weakly to w{x) thus the theorem is proved. □ 

Remark 5.4. If the assumptions in Theorem 4.3 are all satisfied, then it is 
easy to see that as m ^ 00, 7^"*^ converges weakly to the semicircle law 
w{x). 

Remark 5.5. Por a block Hankel matrix as we discussed in Remark 3.3 wc 
also have that as m — )■ 00, 7^^ converges weakly to the semicircle law w{x) 
because (see (3.9)) 

r(m,7r) = tl{ %) | 1 < ^a^, < "i; *ai = ^(6^+1); *(ai+l) = (1 < ^ < 

+ 0(m'=) 

= m2^-^W + 0(m^) 
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and {vr € 7-'2(2fc)|7r is noncrossing} C (2^) and the fact that the integral 
in (3.8) is 1 when vr is noncrossing (see [4]). If the blocks of Hjsi are Hermitian 
matrices, the results would be similar. If the blocks in a block Toeplitz 
matrix T = (^i-j)f^=i are complex matrices and A^g = (^s)"^, then the 
results would also be similar. 

Remark 5.6. If the blocks in a block Toeplitz matrix T = {Ai-j)f-^^ are 
symmetric matrices and A^g = Ag, just like those discussed in [18], then the 
eigenvalue distribution of T jyJmN converges almost surely to a distribution 
T^^"*-* which is determined by its even moments 

7reP2(2/c) ^[0,1]X[-1,1] j=i g=i 

where r{m,7r) is the same as (3.9). So we also have that as m — )• oo, 7^^™'* 
converges weakly to the semicircle law w{x). 
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